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ABSTRACT:  Robust  methods  to  tune  the  unique  electronic 
properties  of  graphene  by  chemical  modification  are  in  great 
demand  due  to  the  potential  of  the  two  dimensional  material  to 
impact  a  range  of  device  applications.  Here  we  show  that  carbon 
and  nitrogen  core-level  resonant  X-ray  spectroscopy  is  a  sensitive 
probe  of  chemical  bonding  and  electronic  structure  of  chemical 
dopants  introduced  in  single-sheet  graphene  films.  In  conjunction 
with  density  functional  theory  based  calculations,  we  are  able  to 
obtain  a  detailed  picture  of  bond  types  and  electronic  structure  in 
graphene  doped  with  nitrogen  at  the  sub-percent  level.  We  show 
that  different  N-bond  types,  including  graphitic,  pyridinic,  and 
nitrilic,  can  exist  in  a  single,  dilutely  N-doped  graphene  sheet.  We  show  that  these  various  bond  types  have  profoundly  different 
effects  on  the  carrier  concentration,  indicating  that  control  over  the  dopant  bond  type  is  a  crucial  requirement  in  advancing 
graphene  electronics. 
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The  two  dimensional  (2D)  nature  of  graphene  endows  it 
with  a  number  of  unique  properties  that  make  it  a 
promising  electronic  contact  material,  including  high  con¬ 
ductivity  (R  ~  20—30  Ohm/square),  transparency,1  and 
tunable  electron  density."'3  The  opportunity  to  tune  the  work 
function  and  carrier  concentration  of  graphene  enables  a  range 
of  electronic  and  optical  applications,  and  diverse  methods  to 
achieve  this  are  being  actively  explored.  Electrostatic  gating  is 
one  well-established  method,"  but  chemical  modification  on 
graphene’s  basal  plane  can  potentially  achieve  the  required 
control  without  the  use  of  external  voltages.  Several  recent 
experiments  have  explored  chemical  modification  of  graphene 
with  foreign  atoms  including  nitrogen,5-14  boron,15'16 
fluorine,  '  and  hydrogen.  In  general,  these  atomic  dopants 
can  be  incorporated  via  several  different  bond  types  in  a 
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graphene  sheet.  A  microscopic  understanding  of  the  different 
bond  types  and  their  diverse  effects  on  the  workfunction,  carrier 
concentration,  and  the  local  electronic  structure  in  nitrogen- 
doped  graphene  films  is  therefore  highly  desirable. 

The  high  energy  resolution  and  tunable  polarization  and 
energy  of  synchrotron  light  enable  the  detection  of  chemically 
distinct  species  and  the  identification  of  bond  types  even  of 
sub-percent  level  dopants  in  monolayer  graphene.  Recent 
scanning  tunneling  microscopy  (STM)  and  spectroscopy 
measurements  have  provided  important  information  on  the 
electronic  structure  of  the  films  near  the  Fermi  energy. 
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However,  STM  is  restricted  to  low  (<1  V)  energy  measure¬ 
ments  so  that  information  on  the  energy  and  symmetry  of  the 
valence  orbitals  and  the  chemical  bond  environments  of  atoms 
is  limited.  Here  we  use  photoelectron  (XPS)  spectroscopy, 
carbon  and  nitrogen  core-level  X-ray  absorption  (XAS)  and 
emission  (XES)  spectroscopy  in  combination  with  density 
functional  theory  (DFT)  based  calculations  to  obtain  an  atom- 
specific  picture  of  dilutely  doped,  single-layer  graphene  (SLG), 
despite  the  small  probability  of  radiative  core  hole  decay  and 
the  weak  signal  from  subpercent  atomic  dopant  levels.  XAS 
provides  direct  information  on  the  dopant  bond  type, 
orientation,  and  concentration.  Coupling  XAS,  which  probes 
the  local  unoccupied  density  of  states  for  a  given  atom  of 
interest,  with  XES,  which  provides  complementary  atom- 
specific  and  symmetry-resolved  information  about  the  local 
occupied  density  of  states,  allows  us  to  identify  the  influence  of 
different  N-dopants  on  the  graphene  electronic  structure. 

Figure  1  shows  angle-dependent  XAS  measured  in  total 
electron  yield  (TEY)  mode  at  the  carbon  and  nitrogen  K-edges 


Figure  1.  C  Is  (left)  and  N  Is  (  right)  XAS  measurements  for  (a) 
pristine  and  (b)  N-doped  graphene  (grown)  on  Cu  foil,  (c)  Strong 
dependence  of  the  integrated  XAS  intensity  (area)  of  the  n*  resonance 
on  incident  angle,  which  is  largest  when  the  E-vector  of  the  light  is  out- 
of-plane  (grazing  angle  10°),  clearly  shows  the  planar  orientation  of 
both  the  C=C  and  C=N  bonds  in  the  graphene  lattice,  (d)  DFT- 
computed  XAS  for  a  substitutional  nitrogen  dopant  bonded  to  three 
carbons  (graphitic  (N3))  in  the  graphene  lattice;  the  match  to 
experiment  is  clear. 


for  pristine  (PC)  and  NH3-grown  graphene  (NGl)  on  Cu  foil 
(sample  preparation  details  in  the  Supporting  Information).13 
The  C  Is  XAS  (left  panel)  shows  the  characteristic  spectral 
signature  of  graphene,  including  strong  peaks  at  ~285  and 
~292  eV,  corresponding,  respectively,  to  the  Is  to  jt*  and  a* 
C=C  bond  resonances,  as  well  as  the  sharp  excitonic  feature  at 
291.5  eV  indicative  of  long-range  order  in  the  electronic 
structure  of  pristine  graphene;  these  spectral  features  are  also 
found  in  graphite.  '  The  spectra  display  the  strong  angular 
dependence  expected  for  the  planar  bonds  of  a  2D  system.19 
The  C  Is  XAS  is  nearly  identical  for  the  pristine  and  N-doped 
graphene  on  Cu  foil,  indicating  that  the  dilute  nitrogen 


concentration  in  the  doped  sample  does  not  disrupt  the  high 
quality  of  the  graphene  sheet.  Importantly,  such  a  disruption 
can  be  observed  for  higher  ammonia  concentrations  during 
CVD  growth,  indicative  of  sample  quality  (Supporting 
Information  Figure  S2). 

In  the  N  Is  XAS  (right  panel),  no  distinct  peaks  appear  in 
the  pristine  graphene  spectra  as  expected.  The  N-doped  sample 
shows  sharp  XAS  peaks  at  ~401  and  408  eV,  corresponding  to 
Is  to  n *  and  <7*  transitions,  respectively,  for  a  single,  well- 
defined  substitution  type.  The  angular  dependence  of  the  N  Is 
XAS  follows  that  of  the  carbon  XAS,  indicating  that  the  features 
are  due  to  planar  N— C  bonds  in  the  graphene  sample  (Figure 
lb,d).  On  the  basis  of  these  two  observations  and  previous 
studies  on  related  systems,  we  assign  the  nitrogen  species  to 
sp2-bonded  graphitic  nitrogen  with  three  C  neighbors."1-"4  For 
comparison,  calculations  based  on  DFT  are  used  to  probe 
spectral  features  for  specific  dopant  structures25,26  (details  of 
the  methods  are  provided  in  the  Supporting  Information 
including  Figure  Si).  The  excellent  match  of  the  computed 
XAS  for  the  N3-dopant  (Figure  Id)  corroborates  this 
assignment,  establishing  a  clear  spectral  fingerprint  of  graphitic 
nitrogen. 

XPS  is  the  determinative  analysis  tool  that  provides 
information  on  atomic  concentration  and  local  coordination 
environment  and  has  been  applied  to  the  study  of  N-doping  in 
graphene.5,9  Figure  2a  shows  the  N  Is  XPS  for  the  same  PG1 
and  NGl  samples  on  Cu  foil  shown  in  Figure  1,  as  well  as  for 
N-doped  graphene  (grown  at  a  lower  NH3-partial  pressure) 
transferred  to  Si02  (NG2/Si02).  NGl  shows  a  weak,  broad 
peak  at  ~400  eV  binding  energy  (BE)  that  we  unequivocally 
assign  to  graphitic  (N3)  dopants  based  on  the  XAS  (Figure  l), 
which  is  in  agreement  with  studies  of  modified  carbon 
films.  '  ’  Using  tabulated  photoelectron  cross  sections  of 
nitrogen  and  carbon,29  a  total  nitrogen  concentration  of  ~0.4% 
is  obtained  for  N-doped  graphene  on  both  the  Cu  foil  and 
transferred  to  Si02.30  In  Figure  2c,  individual  N3  dopants  are 
visualized  for  the  NGl  sample  using  STM,  observed  as  bright 
triangular  spots.  By  simply  counting  bright  spots  in  the  image, 
we  estimate  ~0.3%  substitutional  nitrogen  in  the  graphene 
lattice,  which  is  in  good  agreement  with  the  value  obtained 
from  XPS. 

Whereas  the  nitrogen  dopants  in  NGl  can  be  assigned  to 
graphitic  (N3)  nitrogen,  the  asymmetric  peak  of  the  transferred 
sample  NG2/Si02  (Figure  2a)  indicates  that  (at  least)  two  N— 
C  bond  species  contribute  to  the  total  nitrogen  concentration. 
We  use  Nl,  N2,  and  N3  to  denote  nitrogen  with  one  (nitrile¬ 
like),  two  (pyridinic)  and  three  (graphitic)  carbon  bonds, 
respectively.  In  addition  to  the  N3  peak  (~400  eV),  the 
transferred  NG2/Si02  sample  also  has  intensity  at  lower  BE 
(~398.5  eV),  indicating  less  electronegative  N— C  bonds,  that 
in  previous  studies  have  been  assigned  to  pyridinic  (N2)21’31 
and  nitrilic  (Nl)  groups.21  However,  the  multitude  of  possible 
N-bond  configurations  (graphitic,  pyridinic,  pyrrolic,  nitrilic) 
with  similar  BE  and  the  broadness  of  the  photoemission  peaks 
make  the  identification  of  bond  types  difficult,  as  evidenced  by 
the  range  of  interpretations  of  similar  data  sets."1,31  33 

We  return  to  N  Is  XAS  to  identify  the  origin(s)  of  the  low 
BE  feature  in  the  N  Is  photoemission.  Figure  3  compares 
angle-dependent  N  Is  XAS  measured  in  partial  electron  yield 
(PEY)  mode  for  N-doped  graphene  grown  on  Cu  foil  (NGA, 
Figure  3a)  and  quickly  loaded  into  UHV  with  samples 
transferred  to  Si02  (NGB/Si02,  Figure  3b)  or  on  Cu  foil 
mechanically  bent  during  brief  exposure  to  ambient  atmosphere 
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Figure  2.  (a)  N  Is  and  (b)  C  Is  XPS,  for  pristine  (PGl)  and  N-doped  graphene  (NGl)  samples  grown  on  Cu  foil  and  for  samples  transferred  to 
Si02  (PG2(/Si02),  NG2(/SiOz)).  XAS  for  samples  PGl  and  NGl  are  shown  in  Figure  1.  Quantitative  XPS  analysis  is  based  on  spectra  normalized 
on  the  high  kinetic  energy  side;  the  C  Is  XPS  shown  here  is  normalized  to  maximum  peak  height  to  facilitate  comparison  of  peak  shape,  (c)  STM 
image  showing  individual  graphitic  dopants  in  sample  NGl. 


Figure  3.  Experimental  N  Is  XAS  (gray  curves)  measured  in  PEY 
mode  for  (a)  NGA  (on  Cu);  (b)  NGB  transferred  to  Si02;  (c)  NGC 
(on  Cu)  exposed  to  ambient  atmosphere,  and  DFT-computed  spectra 
(red  curves)  for  cluster  models  of  graphitic  (N3);  pyridinic  (N2);  and 
nitrile-like  (Nl(+H))  nitrogen  dopants,  used  to  assign  the 
experimental  resonances.  Solid  and  dashed  lines  indicate  out-of- 
plane  and  in-plane  intensity,  respectively,  in  the  computed  spectra;  the 
measurements  were  taken  at  grazing  incidence  of  20°  (solid  line)  and 
70°  (dashed  line). 


(NGC,  Figure  3c)  with  the  latter  two  showing  a  low  BE 
component  in  addition  to  the  graphitic  (N3)  peak  (see  Figure 
2a).  Because  the  XPS  BE  is  referenced  to  the  sample  Fermi 
energy,  for  a  metallic  system  such  as  graphene  it  can  be  directly 
correlated  to  the  onset  energy  for  the  XAS  (see,  e.g.,  ref  34). 
Accordingly,  new  n*  resonances  appear  at  lower  photon  energy 
(399  and  399.5  eV)  for  NGB  (on  Si02)  and  NGC  (on  Cu)  in 


addition  to  the  graphitic  peak  (NGA),  each  corresponding  to  a 
distinct  N— C  bond  species  (Figure  3  top).  On  the  basis  of 
DFT-computed  spectra  for  different  N-bond  structures  shown 
in  the  lower  panel  of  Figure  3,  we  assign  the  peaks  at  399  and 
399.5  eV  to  pyridinic  (N2)  and  hydrogenated  Nl-type  dopants, 
respectively,  which  is  in  agreement  with  previous  studies  of 
modified  carbon  systems  (Figure  3  inset).2  ,24'31-33  In  contrast 
to  the  sharp  n*  and  a *  resonances  of  substitutionally  N-doped 
graphene  (NGA),  NGB  and  NGC  show  a  weaker,  broader,  and 
energy-shifted  a*  resonance  as  predicted  by  theory  for  N2  and 
N1  dopants  (Figure  3b,  bottom),  which  is  also  similar  to  the 
case  of  high  NH3  growth  (Supporting  Information  Figure  S2). 
(We  note  that  the  computed  XAS  spectrum  and  electronic 
structure  effect  of  one  pyridinic  (N2)  dopant  at  a  vacancy  site  is 
very  similar  to  the  case  of  three  such  dopants  (see  Table  l);  we 
show  the  case  of  a  single  N2  dopant  here  for  simplicity.) 


Table  1.  Computed  Work  Function,  Doping  Effect,  and 
Fraction  of  Electrons  Added  or  Withdrawn  (— )  from  the 
Lattice  Per  N-Dopant  for  the  Different  Bond  Types 
Compared  to  Pristine  Graphene  (PG)  and  Graphene  with 
Vacancy  Defects 


bond  type 

work  function  (eV) 

doping  effect 

nJN 

PG 

4.43 

- 

- 

vacancy 

4.74 

P 

- 

vacancy  +  H 

4.71 

P 

- 

graphitic  (N3) 

3.98 

n 

0.54 

pyridinic  (N2) 

4.83 

P 

-0.45 

pyridinic  (+H) 

4.29 

n 

0.02 

nitrilic  (Nl) 

4.92 

P 

-0.65 

nitrilic  (+H) 

4.83 

P 

-0.66 

nitrilic  (+2H) 

4.82 

P 

-0.48 

How  do  the  different  N-bond  environments  affect  the 
macroscopic  electronic  structure  of  doped  graphene?  With  clear 
correlations  between  the  spectroscopic  signatures  in  Figures  1 
and  3  and  local  dopant  bonding  motifs,  DFT  calculations  for 
extended  graphene  sheets  with  1%  of  the  dopant  structures35 
(details  provided  in  the  Supporting  Information)  are  used  to 
probe  changes  in  the  work  function  for  the  different  N-bond 
environments  and  the  fractional  charge  added  (withdrawn)  per 
N  atom  for  each  dopant  type  (Table  l).  The  different  N-bond 
types  have  profoundly  different  macroscopic  electronic 
structure  effects  on  the  graphene.  We  find  that  graphitic 
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nitrogen  (N3)  n-dopes  graphene,  consistent  with  the  expected 
picture  of  substitutional  N-doping.13  In  contrast,  pyridinic  (N2) 
and  nitrile-like  (Nl,  with  or  without  H)  dopants,  which  occur 
at  vacancy  sites  and  domain  edges,  and  vacancy  defects 
themselves,  produce  p-type  doping  of  graphene.  Table  1  shows 
that  ~0.5  electron  per  N-dopant  is  added  to  the  carbon  n- 
network  in  the  N3  case,  which  is  in  excellent  agreement  with 
scanning  tunneling  spectroscopy  data.13  On  the  other  hand, 
pyridinic  and  nitrile-like  nitrogen  remove  charge  from  graphene 
resulting  in  p-type  doping  and  an  increase  of  the  work  function, 
which  was  experimentally  confirmed  by  measuring  the 
workfunction  of  transferred  graphene  with  various  degrees  of 
N2  and  Nl  dopants  relative  to  pristine  graphene  (see 
Supporting  Information  Figure  S3).  This  understanding, 
coupled  with  increasing  control  over  N2/N3  ratios  in  N- 
doped  graphene,14  opens  new  avenues  for  tailoring  the  carrier 
density  and  electronic  properties  of  graphene  at  the  atomic 
level. 

In  the  following,  we  show  that  core-level  spectroscopy  and 
DFT  are  sensitive  to  the  details  of  the  effect  of  the  different 
dopant  bond  structures  and  vacancy  defects  on  the  electronic 
structure  of  graphene.  The  diverse  effects  on  the  carrier  density 
in  graphene  are  visualized  in  computed  charge  density 
difference  (CDD)  plots  in  Figure  4a,  which  show  the 
redistribution  of  charge  induced  by  the  different  bond 
structures  (details  provided  in  the  Supporting  Information). 
Close  to  the  dopant  center,  there  are  substantial  local  changes 
in  the  charge,  including  bond  polarization  effects.  Further  from 
the  dopant  center,  an  increase  in  charge  density  of  the  carbon 
^■-network  is  uniquely  observed  for  the  N3  dopant;  Nl  and  N2 
groups  and  vacancy  defects  on  the  other  hand  withdraw  charge 
from  the  ^-orbitals  (see  Supporting  Information  Figure  S5). 

The  different  N-bond  environments  not  only  result  in  a  clear 
signature  in  the  N  Is  XAS  spectra,  but  at  the  present  doping 
levels,  they  can  also  clearly  influence  the  C  Is  XAS  and  XES 
spectra  (Figure  4b).  Figure  4b  compares  C  Is  XES  and  XAS  for 
pristine  (PG2/Si02)  and  N-doped  (NG2/Si02)  graphene 
transferred  to  Si02  (XPS  shown  in  Figure  2a)  in  which  Nl,  N2 
(p-type),  and  N3  (n-type)  dopants  coexist  (see  Figure  3  and 
Supporting  Information  Figure  S3).  Estimates  of  the 
concentration  of  Nl,  N2,  and  N3  dopants  in  NG2/Si02 
from  XPS  are  provided  in  Supporting  Information  Figure  S4. 
Angle-dependent  C  Is  XAS  shows  that  the  transferred  samples 
are  significantly  rippled  compared  to  graphene  grown  directly 
on  Cu  foil  (see  Supporting  Information  Figure  S5).  The 
deviation  from  planarity,  an  expected  consequence  of  litho¬ 
graphic  transfer  methods,36  results  in  a  mixing  of  a  and  n 
signals  in  the  XES  emission  channels  and  polarization- 
dependent  XAS.  The  in-plane  XAS  and  normal  emission  XES 
are  nonetheless  dominated  by  c-states,  while  the  out-of-plane 
XAS  (and  corresponding  XES)  are  dominated  by  ^-states.  We 
also  note  the  presence  of  small  concentrations  (~0.1%)  of  Nl 
and  N2  dopant  structures  observed  in  PG2/Si02  with  XPS  (see 
Supporting  Information  Figure  S9).  This  will  act  to  reduce  the 
spectral  differences  observed  between  samples  PG2/Si02  and 
NG2/Si02  induced  by  Nl  and  N2  type  dopants,  nonetheless 
Nl  and  N2  bond  signatures  in  sample  NG2/Si02  are  still 
evident. 

Distinct  features  in  the  C  Is  XES  and  XAS,  labeled  1—6  in 
Figure  4b,  are  observed  for  NG2/Si02  relative  to  PG2/Si02. 
These  features  reflect  the  interplay  of  dopant-induced  changes 
to  the  electronic  structure,  which  are  represented  schematically 
in  Figure  4c  based  on  computed  pDOS  and  simulated  XES  for 


Figure  4.  (a)  Charge  density  difference  (CDD)  plots  for  a  single  N3 
(left),  N2  (center),  and  Nl-type  (right)  dopant.  The  plots  show  the 
gain  (red)  and  loss  (blue)  of  charge  in  the  C  ^-network  for  the  n- 
(N3)  and  p-type  (N2  and  Nl)  dopants,  respectively;  arrows  are  drawn 
to  guide  the  eye.  The  CDD  is  taken  relative  to  pristine  graphene  (PG) 
in  the  N3  case  (pdiff  =  p(NG)  —  p(PG))  and  relative  to  PG  with 
vacancy  (CDD  shown  in  Supporting  Information  Figure  S5)  defect(s) 
in  the  N2  and  Nl  cases  (Pag  =  p(NG)  —  p(PGvac))  where  PG^ 
indicates  pristine  graphene  with  a  vacancy,  (b)  C  Is  symmetry- 
resolved  XES  (left)  and  in-  (IP)  and  out-of-plane  (OOP)  XAS  (right) 
for  pristine  (PG2/Si02)  and  N-doped  graphene  (NG2/Si02) 
transferred  to  Si02  for  which  XPS  is  shown  in  Figure  2.  Note  that 
the  buckling  on  the  Si02  generates  a  significant  Ji  and  a  mixing  in  both 
XAS  and  XES;  total  (bottom),  cr-dominated  (center),  and  n- 
dominated  (top)  spectra  are  shown  (see  text).  The  XAS  data  is 
normalized  by  area  past  the  continuum  (330  eV).  (c)  Schematic  of  the 
tr-band  structure  of  PG  (black  curve,  blue  fill  for  occupied  states,  white 
for  unoccupied)  and  distinct  C  Is  core  and  Fermi  level  shifts  and  new 
states  (red  and  green  fill  correspond  to  the  sum  of  occupied  n-  and  a- 
states,  respectively,  for  the  carbon  atoms  in  the  vicinity  of  the  N- 
dopant)  based  on  computed  pDOS  and  simulated  XES  provided  in  the 
Supporting  Information.  C=N,  C=C,  and  CV=C  signify  the  carbon 
bonded  to  the  N-dopant,  the  next-nearest  neighbor,  and  to  the  C  atom 
at  a  vacancy  site  (N2  and  Nl),  respectively.  Specific  values  for  the 
shifts  are  listed  in  the  Supporting  Information. 


graphene  doped  at  the  1%  level  (Supporting  Information 
Figures  S7  and  S8).  First,  there  is  an  overall  shift  of  the  Fermi 
level  with  doping  (see  Table  l).  Second,  the  C  atoms  that  are 
near  the  point  defects  have  a  shift  in  the  local  Is  level  and  a 
change  in  the  local  density  of  states  due  to  interactions  with  the 
N  centered  region.  Since  XAS  and  XES  spectra  result  from 
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atom-localized  absorption  or  emission  events,  these  sites  make 
distinct  contributions  to  the  spectra  that  are  in  aggregate  clearly 
detectable. 

The  N-induced  changes  in  the  ^-dominated  XES/XAS  data 
(Figure  4b  top)  are  well  represented  by  band  structure 
differences  between  N3-doped  graphene  relative  to  PG  (Figure 
4c),  suggesting  that  graphitic  dopants  dominate  this  channel. 
Specifically,  we  observe  a  suppression  of  ^-states,  particularly 
around  the  peak  in  the  density  of  states  below  the  Fermi  level 
Ef  resulting  in  the  suppression  of  intensity  in  the  280—282  eV 
range  (feature  4  in  Figure  4b  top).  The  calculations  also 
indicate  a  shift  of  the  center  of  ^-states  near  the  N3  dopant  to 
higher  BE  (lower  photon  energy),  and  development  of  a 
significant  peak  in  the  local  density  of  states  at  the  bottom  of 
the  n  band  (~9  eV  below  EF,  see  Supporting  Information 
Figure  S7).  These  changes  induced  by  the  N3  dopant  lead  to 
the  increase  of  intensity  in  the  XES  below  278  eV  (feature  6  in 
Figure  4b  top).  Moreover,  a  shoulder  appears  on  the  low 
energy  side  of  the  n*  resonance  at  285.5  eV  of  the  NG2/Si02 
XAS  (feature  1  in  Figure  4b  top);  we  note  that  a  corresponding 
shoulder  appears  on  the  low  BE  side  of  the  C  Is  photoemission 
for  this  sample  (Figure  2b).  This  shoulder  results  primarily 
from  partially  unoccupied  carbon  ^-states  due  to  N3,  see  Figure 
4c.  This  is  consistent  with  a  recent  report  by  Schultz  et  al.3  in 
which  the  pre-edge  shoulder  on  the  C=C  n*  resonance  was 
attributed  to  doping,  rather  than  to  vacancy  or  edge  defects37 
and  shown  to  be  independent  of  the  degree  of  corrugation  in 
the  transferred  layer,  which  shifts  the  graphene  Dirac  point.  We 
note  that  carbon  atoms  at  vacancy  defects  (with  or  without  N 
present)  may  contribute  to  this  shoulder,  although  to  a  much 
lesser  degree  (see,  e.g.,  N1  in  Figure  4c  and  Supporting 
Information  Figure  S7).  (The  simulated  spectra  for  N1  and  N2 
dopants  show  a  similar  suppression,  as  well  as  a  modest 
enhancement  in  the  282—284  eV  range  (Supporting 
Information  Figure  S8)  that  is  not  observed  in  the  high  energy 
tail  of  the  XES.) 

On  the  other  hand,  the  changes  in  the  normal  emission  (<7- 
dominated)  XES  for  NG2/Si02  relative  to  PG2/Si02  are 
largely  represented  by  the  electronic  effects  of  the  N2  and  N1 
dopants.  All  three  dopants  show  a  suppression  of  a  local 
density  of  states  in  the  range  5—7  eV  below  the  Fermi  energy 
and  an  overall  redistribution  of  the  density  of  states  to  lower 
energy  (Supporting  Information  Figure  S7),  which  contributes 
to  the  reduced  intensity  around  278  eV  (feature  5  in  Figure  4b 
middle).  However,  N2  and  N1  dopants  introduce  states  with  a 
character  close  to  the  Fermi  energy  (Figure  4c  and  Supporting 
Information  Figure  S7),  which  correspond  to  features  2  and  3 
in  Figure  4b,  center;  the  simulated  XES  clearly  shows  added 
intensity  around  282  and  280  eV  for  N2  and  N1  respectively  in 
the  normal  emission  signal  (Supporting  Information  Figure 
S8).  N2  and  N1  (without  hydrogen)  introduce  occupied 
nonbonding  (lone-pair)  (7-states,  while  the  hydrogenated  N1 
dopants  (Figure  3)  introduce  occupied  N— H  (7-states  (see 
pDOS  in  Supporting  Information  Figure  S7).  Interestingly, 
simulations  show  that  the  isolated  vacancy  introduces  addi¬ 
tional  intensity  near  283  eV  (Supporting  Information  Figure 
S8),  which  may  also  contribute  to  the  significant  enhancement 
of  intensity  observed  in  the  high  energy  tail  of  the  normal 
emission  XES  (Figure  4b  middle). 

In  addition  to  distinct  EF  shifts  for  the  different  N-dopant 
bond  structures  (Table  l),  the  C  Is  level  of  the  carbon  atom 
bonded  to  the  N-dopant  and  that  of  the  next-nearest  neighbor 
(C=C)  are  also  sensitive  to  the  details  of  the  N-bond 


environment  and  significant  variation  in  the  magnitude  and 
direction  of  the  shift  is  observed  across  the  N-dopant  bond 
series  (values  provided  in  Supporting  Information).  The  C  Is 
level  of  the  carbon  atom  at  a  vacancy  (CT=C)  shifts  (by 
different  amounts)  to  higher  BE  for  the  N2  and  Nl-type 
dopants  (Figure  4c);  interestingly,  the  opposite  shift  is 
observed  for  a  carbon  atom  at  vacancy  site  in  pristine  graphene. 
(The  calculations  indicate  that  the  electronic  structure  of 
graphene  is  only  perturbed  in  the  close  proximity  (first  two 
coordination  shells)  of  the  N-dopant,  as  visualized  in  Figure  4a 
and  consistent  with  our  previous  findings  for  individual  N3- 
dopants.13)  The  broader  energy  distribution  of  the  occupied  n- 
and  (7-states  (~0— 6  eV  relative  to  EF),  as  well  as  the  distinct 
core-level  and  Ep  shifts  for  Nl,  N2,  and  N3  dopants  (Table  1 
and  Figure  4c),  account  for  the  width  of  the  shoulder  in  the 
measured  XES,  while  the  partially  unoccupied  ^-states  produce 
the  XAS  shoulder. 

The  observed  electronic  effect  of  the  different  bond 
structures  in  N-doped  graphene  can  be  understood  in  a  simple 
chemical  bond  picture,  summarized  in  the  following.  In  the  N3- 
case,  four  of  the  five  valence  electrons  in  N  fill  the  a-  and  n- 
orbitals  as  for  C,  leaving  one  extra  electron.  The  additional 
charge  is  to  ~50%  localized  on  the  N3  dopant,  which 
electronically  couples  to  its  nearest  C  neighbors.  Consequently 
we  observe  a  single  N  nonbondin^  ^-orbital  resonance  in  the 
projected  density  of  states  (pDOS)  3  and  in  the  simulated  XES 
(Supporting  Information)  at  the  Fermi  energy,  and  a  largely 
unoccupied  pz  state  on  the  N3-bonded  carbon  (Figure  4c).  The 
remaining  charge  of  ~0.5  electron  per  N3  dopant  is  distributed 
in  the  local  network  of  carbon  ^-states  and  n-dopes  the 
system;13  this  appears  as  a  gain  in  orbital  occupation  (pDOS) 
and  charge  density  of  the  7T-orbitals  (Figure  4a  left).  In  contrast, 
N2  and  Nl  dopants  have  the  opposite  electronic  effect;  they 
withdraw  charge  from  the  carbon  ^-network  (Table  1  and 
Figure  4a)  to  p-dope  graphene.  In  the  N2-case,  two  electrons 
fill  ct-bonds  with  carbon  neighbors,  two  electrons  form  a  lone 
pair  in  the  graphene  plane,  and  the  remaining  electron  occupies 
the  N  estate  (Figure  4c  and  Supporting  Information).  Thus, 
the  N2  dopant  has  the  equivalent  occupation  of  a  nominal 
carbon  in  graphene,  but  a  ^-electron  is  missing  due  to  the 
vacancy  site  so  the  system  is  p-doped.  Similarly,  the  Nl  dopant 
forms  one  tr-bond  with  its  C  neighbor,  forms  (7-bonds  to  two  H 
(or  one  N-lone-pair  orbital  in  the  unhydrogenated  case).  Two 
electrons  go  into  the  ^-network,  but  as  for  the  N2  case,  the 
missing  ^-electrons  from  the  divacancy  sites  that  accommodate 
this  form  of  the  N  dopant  result  in  overall  p-type  doping,  which 
is  further  enhanced  by  the  polarization  of  charge  in  the  plane  of 
the  N-dopant  (see  Figure  4a  center).  In  contrast,  if  the 
conditions  for  the  formation  of  the  lone-pairs  or  other  localized 
states  (like  edge  states)  are  unfavorable,  for  example,  in  the 
graphitic  case  or  hydrogenated  pyridinic  groups  (Table  l),  the 
additional  charge  from  nitrogen  is  forced  to  go  to  the  extended 
carbon  ^-network,  essentially  neutralizing  the  p-doping  effect  in 
the  N2  case  (Table  l). 

In  summary,  we  have  presented  a  microscopic  understanding 
of  the  different  bond  types  and  their  diverse  effects  on  the 
workfunction,  carrier  concentration  and  the  local  electronic 
structure  in  nitrogen-doped  graphene  films  from  core-level 
spectroscopy  and  DFT.  These  findings  show  that  controlling 
the  bonding-type  of  dopant  substitutions  in  graphene  is 
essential  when  designing  properties  of  graphene  as  functional 
component  in  future  devices.  This  work  demonstrates  that 
synchrotron  radiation  based  spectroscopies  provide  invaluable, 
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atom-specific  tools  to  determine  the  resulting  electronic 
properties  of  different  dopant  and  defect  structures  in  graphene, 
even  at  subpercent  dopant  levels.  Combined  with  continuing 
advances  in  the  growth,  controlled  doping  of  graphene,  and 
cleaner  transfer  procedures,  the  detailed  correlation  of  bond 
type  and  electronic  structure  demonstrated  here  promises  to 
facilitate  the  atomic-level  control  of  electronic  properties 
desired  for  next  generation  graphene-based  devices. 
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